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<N ■ ABSTRACT 

We present spectropolarimetric radio images of the supernova remnant (SNR) 
. G296. 5+10.0 at frequencies near 1.4 GHz, observed with the Australia Telescope 

Compact Array. By applying rotation measure (RM) synthesis to the data, a 
pixel- by-pixel map of Faraday rotation has been produced for the entire remnant. 
We find G296. 5+10.0 to have a highly ordered RM structure, with mainly positive 
r£3 ■ RMs (mean RM of +28 rad m -2 ) on the eastern side and negative RMs (mean RM 

of —14 rad m~ 2 ) on the western side, indicating a magnetic field which is directed 
away from us on one side and toward us on the other. We consider several possible 
mechanisms for creating the observed RM pattern. Neither Faraday rotation 
in foreground interstellar gas nor in a homogeneous ambient medium swept up 
by the SNR shell can easily explain the magnitude and sign of the observed 
RM pattern. Instead, we propose that the observed RMs are the imprint of an 
azimuthal magnetic field in the stellar wind of the progenitor star. Specifically, we 
^ ! calculate that a swept-up magnetized wind from a red supergiant can produce 



6 
N 



RMs of the observed magnitude, while the azimuthal pattern of the magnetic 
field at large distances from the star naturally produces the anti-symmetric RM 
pattern observed. Expansion into such a wind can possibly also account for the 
striking bilateral symmetry of the SNR's radio and X-ray morphologies. 

Subject headings: ISM: individual (G296. 5+10.0), magnetic fields — polariza- 
tion — stars: winds — supernova remnants — radio continuum: ISM 
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Introduction 



Observations of radio supernova remnants (SNRs) provide a unique probe of the cir- 
cumstellar environments of supernova progenitor stars and the interstellar medium (ISM). In 
particular, the interactions between SNRs and their environments give rise to a wide range of 
observable effects from which we can derive important physical insights. Spectral line maps 
of neutral hydrogen (Hi) in the vicinity of SNR s enable us to calcu l ate the age, swept-up 



shell mass and ene rgy of the supernova explosion ( jGaensler et al.lll998t IVelazquez et al.ll2002 



Foster et al.l I2004J ) . Combined with model rotation cur ves for the Galaxy, Hi data may be 
used to constrain the kinematic distance to the remnant (IPineault et al.lll993t iTian fc Leahy 
20081 ). Under certain conditions, regions of interaction between a SNR shock and a dense 
molecular cloud produce hydroxyl masers, allowing us to derive the properties of the SNR 
shock, estimate the distanc e to the remnant and to calculate the postshock magnetic field 
through the Zeeman effect ( IFrail et al.lll994j ; IClaussen et al.l 119971 ). In many SNRs there is 
evidence of physical confinement by denser, cooler surrounding material, w hich allows us 



to probe the physical, chemical and magnetic properties of molecular clouds (IBurton et al. 



19881 ; IFrail et al.lll994j ; iGreen et al.lll997l ; IChevarierlll999l) and provide information on parti 



cle acceleration and the production of cosmic rays (e.g.. iKoyama et al.lll995t iTanimori et al. 
19981 : iFukui et al.ll2003h . 



The appearance of young SNRs can be strongly influenced by the mas s- loss history of 
the progenitor star ( jChevalierl Il982l ; iFranco et al.lll991t iKoo fc Heilesl Il995h . For example, 
the properties of the progenitor stellar wind have been inferred in the young SNR Cas- 
siopeia A by studying the dynamic s of the clumpy shock boundary in the remnant and 
calculating the swept-up shell mass ([Chevalier &: Oishill2003t IVink et al.lll998l ). In addition, 
the structure and brightness evolution of the radio remnant of supernova 1987A has been 
shown to have been shaped by successive interactions of the remnant with material from 
stellar winds originating from the blue- and red -supergiant phases of the progenitor star 
( jChevalier fc Franssonlll987l ; iGaensler et al.l 119971 ). Interaction of the supernova ejecta with 
stellar winds is a phenomenon that occurs in the earlier or mid-stages of a radio SNR's life 
( )Dwarkadasll2005l ) although in some cases the passage of the remnant through the circumstel- 
l ar material (CSM) can continue for many thousands of years in large stellar wind bubbles 
( ILandecker et al.l Il999l ; iGvaramadzd 120061 ) . 



Magnetic fields a re known to influence the development of supernova remnants ( IKulsrud et al 



1965c IChevalierlll974j ) . This can be as a direct consequence of magnetic c onfinement of the ex 
panding SNR shell under th e influence of external ISM magnetic fields ( Ivan der Laanlll962a 
Whiteoak fc Gardner! Il968l ). or due to the sweeping-up of a CSM whose distribution is al 
ready anisotropic due to magnetic channeling and confinement of the pre-supernova wind(s 
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see, e.g., lud-Doula fc Owockil 12002 : lud-Doula et al.ll2008f l 



Alter natives to the direct shaping o f SNRs by magnetic confinement have also been 
proposed. iBisnovatyi-Kogan et al.l ( 119901 ) asserted that the eventual shape of a supernova 
remnant is primarily influenced by the shape of the cavity swept out by the progenitor star. 
In some cases, the shape of such a stellar w ind bubble will b e domi nated by magnetic fields. 
In a theoretical treatment of this problem, IChevalier fc Lud (119941 ) considered the case of a 
bubble driven by a magnetized stellar wind. They found that, in a bubble resulting from a 
stellar wind with a toroidal magnetic field, the magnetic tension in the equatorial region of 
the shell can confine the bubble and cause elongation in the polar direction. The proposal 
that S NRs can be shaped by a magnetically-influenced cavity was arrived at by iGaensler 
(119981 ). who proposed that SNRs with barrel morphologies (remnants with negligible emission 
along the axis of symmetry and bright limbs of emission perpendicular to this axis) form in 
ISM bubbles that are elongated by the Galactic magnetic field. Whether this elongation is 
due directly to magnetic confinement or the product of density stratifications along the local 
magnetic field axis is less clear. 

All of these different mechanisms will have an effect at some point in the life of a SNR. It 
is likely that very young remnants are dominated by the supernova explosion itself, middle- 
aged remnants will primarily be sculpted by the CSM and older remnants will be shaped by 
the surrounding ISM. Disentangling these different effects is a complicated task, but if the 
production and evolution of supernova remnants is to be understood, quantitative studies of 
the magnetic fields in their immediate environments are crucial. 

A useful experimental probe of cosmic magnetic fields is Faraday rotation. When polar- 
ized electromagnetic radiation propagates through a magnetized ionized medium, its electric 
field vector is rotated through an angle which depends on the square of the wavelength 
and the path integral of the electron- weighted line-of-sight magnetic field strength. This 
effect is called Faraday rotation. The angle of rotation, A%, is related to the square of the 
wavelength, A, via A% = RMX 2 where the RM is the rotation measure and is defined as 



RM = 0.81 



observer 



n P B\\ dl rad m 2 



where B\\ is the line-of-sight component of the magnetic field strength in /xG, n e is the 
electron density in cm -3 and dl is measured in parsecs. 

Evaluating the RM of linearly polarized radio emission gives us an indication of the 
line-of-sight magnetic field orientation in the foreground Faraday rotating medium. Mea- 
surements of other physical cues, such as the interstellar rotation measure along similar lines 
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of sight, can assist us in disentangling the Faraday rotation due to the Galactic foreground 
from the RM occuring internally to the source being studied. If the RM is intrinsic to the 
source, by measuring RMs at many points in an image we can learn about the line-of-sight 
orientation of the smooth component of the magnetic field in that source. In a supernova 
remnant, this will either tell us about the ambient magnetic field in the ISM, or the magnetic 
field in the swept-up shell of circumstellar material depending on the age and history of the 
remnant. 

In this paper we consider the Faraday rotation against G296. 5+10.0 (also known as 

PKS 1209-51/52), a high-Galactic latitude supernova remnant with a barrel morphology at 

both radio and X-ray wavelengths (Whiteoak & Gardner 1968; Roger et al. 1988; Helfand & Becker 
19841 ) . The limbs of the remnant are believed to mark the material swept up by the blast 
wave from a core-collapse supernova as evidenced by the radio-quiet X-ray emitting neu 
tron star (IE 1207.4-5209) at its center jHelfand k Becker! Il984 



G296. 5+10.0 displays a high degree of reflection symmetry ( jStorey et al.lll992|). and lies 1 



Zavlin et al. 



degrees above the Galactic plane at an estimated distance of 2.1 kpc (IGiacani et al.ll2000l ). 



2000 



2004J). 



Attempts to m easure the polarization and r otation measure i n G2 96. 5+ 10.0 have previ- 
ously been made by I Whiteoak fc Gardner! (119681 ) ; iDickel &; Milnd (119761 ) and lMilne &; Haynes 
(119941 ). but these were limited in their accuracy because the RMs were derived by a linear 
fit to only two or three points in a plot of position angle vs. A 2 . This method of determining 
RMs introduces significant uncertainties into the quoted values, because the polarization 
position angle is k nown only as a multiplicati ve factor of nir radians. By using rotation 
measure synthesis (IBrentjens fc de Bruynll2005l ) to produce a well-sampled rotation measure 
spectrum at each image pixel, it is possible to eliminate these ambiguities and determine 
accurate RMs for the entire region. Another advantage of RM synthesis is that it allows us 
to resolve multiple RM components within a single line of sight. With the flexible correlator 
on the Australia Telescope Compact Array (ATCA) this can be achieved within a single 
128 MHz frequency band. 

In this paper, we present full polarization and rotation measure synthesis images of 
G296. 5+10.0 using RM synthesis of ATCA data. In §2] we introduce the observations and 
methods of data analysis and in §3]we show total intensity, polarization and rotation measure 
images and plots of RM distribution in G296. 5+10.0. The RM distribution and magnetic 
fields in G296. 5+10.0 are discussed in fjU $5] describes our conclusions and suggestions for 
future work in this field. 
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2. Observations & Analysis 



The supernova remnant G296. 5+10.0 was observed between the 8th and the 11th of 
October 1998 using the ATCA . Analysis of the Hi data from these observations was pub- 
lished by iGiacani et al.l (120001 ) . We obtained the data from the ATCA archive in order to 
extract and analyze the (as yet unpublished) commensally recorded multi-channel contin- 
uum polarization information. The observations were carried out in the 210 m configuration, 
resulting in an angular resolution of 1(5 at the central frequency of 1384 MHz. The total 
bandwidth of 128 MHz was correlated in 32 x 4 MHz overlapping spectral channels. Those 
channels affected by self-interference, edge channels and overlapping channels were discarded 
to leave a useable bandwidth of 104 MHz split into 13 x 8 MHz channels. Calibration terms 
were determined using observations of the sources PKS 1215-457 (gain solutions) and PKS 
1934-638 (flux, p olarization leakage an d bandpass solutions) using the software reduction 
package miriad (ISault &: Killeenl 12004 ) . 109 individual pointings were combined to produce 
mosaiced images of the entire remnant in Stokes-/, Q, U, V and polarized intensity. Images 
were deconvolved using the maximum entropy method using the miriad task pmosmem 
( ISault et al.lll999l ). The RMS noise levels in the final continuum images were 2.0, 1.1 and 
1.2 mJy beam -1 for Stokes-/, Q and U respectively. In order to make RM synthesis images, 
we also made separate mosaic images of each individual spectral channel for Stokes- Q and U, 
which had mean RMS noise levels of 3.4 and 3.5 mJy beam -1 , respectively. The maximum 
angular scale to which an interferometer is sensitive is set by the shortest spacing between 
antennas, in this case 31 metres. As a result, any emission greater than 24' is not detected 
in our ATCA observations, resulting in a Stokes-/ flux on large scales below the true value. 
Polarization images of the region are likely to contain almost all the flux however, as differ- 
ential Faraday rotation in the Galactic foreground is expected to break the polarized image 
into smaller-scale structures. 



Rotation measure synthesis (jBrentjens fc de Bruynl 120051 ) was employed to produce a 
full RM spectrum for each image pixel in steps of 20 rad m -2 between —2000 and +2000 
rad m -2 . The RM function had a FWHM of 530 rad m -2 . We then used the spectral 



deconvolution algorithm rmclean ( IHeald et al.l 120091 ) with a maximum of 1000 iterations 
and a clean cutoff of 0.0001 to remove the sidelobe pattern resulting from incomplete A 2 
sampling. The outcome of these procedures was a spectrum of polarized flux as a function 
of RM at each pixel, with typical signal-to-noise ratio of ~30 and a median RM uncertainty 
of 5.5 rad m -2 . Finally, an image of polarized intensity across the remnant was generated by 
using a 3-point polynomial fit around the peak of the deconvolved rotation measure spectrum 
at each pixel. Although we are able to confidently detect a polarized signal at a signal-to- 
noise rat io of ~ 4 — 5, a meaningf ul RM can only be extracted for a somewhat stronger 



signal (jBrentjens fc de Bruynl 120051 ). Accordingly, regions of the final polarized image with 
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a signal below 10 x the root mean squared noise level in that image were masked (excluded) 
from the final maps of RM and polarization. At these high thresholds, polarization bias is 
negligible and so no debiasing was applied to any of the data. 

3. Results 

Figured] presents the images of (a) total intensity, (b) total polarized intensity (the peak 
P as a function of RM in the RM synthesis spectrum) and (c) & (d) peak rotation measure 
for SNR G296. 5+10.0. We see a clear bilateral morphology in both the total intensity and 
the polarized intensity images, with bright limbs of emission to the east and west of the SNR 
symmetry axis. Along this axis, there appears to be a complete absence of radio emission. 
The mean fractional polarization of G296. 5+10.0, found by comparing the flux density of 
the total polarized intensity and the Stokes-/ images at each pixel, is approximately 35%. 
This is likely an upper limit, as our observations are not sensitive to large-scale total power 
emission, as discussed in §2J Figure [T] (c) shows the rotation measures as boxes, denoting 
the RM value averaged over five adjacent pixels. Filled boxes, representing positive RMs, 
dominate the eastern limb and unfilled boxes (negative RMs) domiate the western limb of 
the remnant. Figure [1] (d) shows the same information as an image. Figure [2] shows a 
typical rotation measure spectrum for a pixel in each of the eastern and western limbs of 
G296. 5+10.0. The polarized flux density is similar in each limb of the SNR, but the RM 
changes from positive to negative on either side of the symmetry axis. 

In order to quantify the distributions of rotation measures in G296. 5+10.0 in Figure [3] 
we plot RM against Right Ascension (upper panel) and against Declination (lower panel). 
In both plots, the mean RM in each limb is displayed as a dashed line. The upper panel 
of Figure [3] shows a jump in RM distribution between the limbs of the SNR. The western 
limb has a mean RM of — 14 ± 80 rad m -2 and the eastern limb has a mean of +28 ± 65 
rad m -2 . In contrast, in the lower panel of Figure [3] we see no such gradient or discontinuity 
in RM as a function of declination. This observation is crucial to our understanding of the 
magnetic field in the region, which we discuss further in § §4.2[ 14.31 and 14.41 

4. Discussion 

4.1. Comparison with previous results 

The distribution of RM in Figured] along the limbs of G296. 5+10.0 has a clear symmetry: 
the eastern limb shows predominantly positive RMs, while the western limb displays negative 
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RMs. Since the sign of the RM corresponds to the direction of the electron-weighted average 
magnetic field between the observer and the source of emission, our results provide a high- 
angular resolution view of the orientation of the magnetic field responsible for the Faraday 
rotation. Specifically, the average magnetic field vector points generally toward from us on 
the eastern side, and away from us on the western side. 

We first compare our r esults to previous stud i es of polarization and Faraday rotation 
toward SNR G296. 5+10.0. IWhiteoak fc Gardner! (119681 ) first studied the magnetic fields 
in G296. 5+10.0 using the Parkes 64-m radio telescope to map polarized radio continuum 
emission from the remnant at 0.6, 1.4 and 2.7 GHz. They made maps of RM across the 
source by measuring the best-fit linear relationship between the polarization angle and the 
square of the wavelength at each image pixel. The resulting RM map showed a positive 
range of RMs (+17 to +36 rad m -2 ) in the eastern limb with a positive gradient of RM with 
increasing declination. In the western limb, t hey found a spread o f negative (—4 to —14 
rad m~ 2 ) RMs with no declination dependence. iDickel fc Milnd (119761 ) used a similar method 
to compare the polarized radio continuum emission at 2.7 and 5.0 GHz. They noted that the 
remnant has a moderate east-west RM gradient and that the projected magnetic field has an 
orientation that runs tangential to the bright edges of the SNR. They also measured an east- 
west gradient in rotation measure of +38 to —18 rad m~ 2 w ith uncertainties of ±23 rad m -2 , 



which agree with the values found by IWhiteoak fc Gardner! (I1968I ). iMilne fc Haynesl (119941 ) 
observed G296. 5+10.0 using the Parkes 64-m telescope at 2.4, 4.8 and 8.4 GHz. This study 
surpassed the two previous polarization observations of the region, having both a higher 
frequency and angular resolution than measured previously. They too found strong evidence 
for a tangential plane of sky magnetic field. Their RM maps indicated largely negative RMs 
across the majority of the remnant, but with a large patch of positive RM values around the 
middle of the eastern limb. 



Our re sults agree broadly wit h the east-west gradient of RMs found by lWhiteoak fc Gardner 



(119681 ) and lDickel fc Milnd ( 119761 ). although we find no evidence for a declination dependence 
of RM in the easter n limb . However, our RMs are not readily compared with those derived 
by IMilne &: Haynesl (119941 ); they quote a much higher range of RMs. This may be in some 
part due to their cutoff level for rotation measures being 5% of the peak polarized flux at 
that frequency, which is a much less stringent cutoff criterion than our 10a in the peak 
polarized intensity. The upper limit on percenta ge polarization we deri ve (35%) broadly 
agrees with that ("in excess of 30%") quoted by IMilne fc Haynesl (11994 ). Our data were 
taken around 1.4 GHz, a low radio frequency where extrapolation of measured polarization 
angles to intrinsic values would give large (> 1 radian) uncertainties. We do not provide such 
extrapolations here and assume the tangential distribution of intrinsic projected magnetic 
field vectors calculated by previous authors. 
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Before we interpret the polarization information gained from our observations, we need 
to consider whether the observe d Faraday rotation is predominantly occurring in unrelated 



foreground interstellar gas (e.g.. ICaswell et al.ll2004f) . or in materia l that has been swept up 



by the expanding SNR (IMatsui et al 



1984; 



Kothes fc Brownl 120091 ) . 



4.2. The RM of the Foreground ISM 

Radio pulsars are typically strongl y linearly polarized, an d consequently many such 
sources in the ATNF Pulsar Catalogue! ( Man Chester et al.l 120051 ) have measured RMs. Pul- 



sars are at a variety of distances within the Galaxy, and their Faraday rotation is thought 
to have no intrinsic component due to the pulsar itself. We can thus compare the RMs 
seen against G296. 5+10.0 to the RMs of pulsars in a similar direction and at comparable 
distances, to assess whether the RMs seen toward G296. 5+10.0 are due to the foreground 
ISM. 

There are two pulsars in the catalogue which sit within 10° of G296. 5+10.0, and have 
estimated distances between 1 and 3 kpc (bracketing the distance to the SNR): P SR B1133 



55 U = 292?3, b = +5?9) has RM = +28 ± 5 rad m~ 2 at a distance of » 2.6 kpc flQiao et al. 



19951 ) while PSR J1253-5820 (£ = 30 3?2, b = +4?5) has RM = +18 ± 1 rad m~ 2 at a 



distance ~ 2.9 kpc ( Noutsos et all2008h. whe r e the distance estimates come from the NE2001 



Galactic electron density of ICordes fc Laziol (120021 ). If we assume that the observed RM is 



the projection of a magnetic field predominantly parallel to the Galactic plane, and that the 
field is sufficiently uniform along these sightlines that the magnitude of the RM scales with 
distance, we can estimate from this pulsar data that the RM of the foreground ISM between 
Earth and the SNR is approximately +10 to +20 rad m -2 . 

We do not find a good match between this expectation and the observed distribution of 
RMs against the SNR. The RMs of the western limb are of the opposite sign to the pulsar 
RMs, while those in the east are of the correct sign, but typically have a magnitude twice as 
large as what would be expected from the pulsar data. 

Additional evidence against a foreground interpretation for the SNR RMs is the relative 
uniformity of the sign of the RM distribution along each limb (shown in the lower panel of 
Figure [3]), but the difference in the sign of RM between limbs (upper panel of Figure E]). If 
this Faraday rotation occurs in the foreground, it requires a complete reversal of the mag- 
netic field on a relatively small scale in the east-west direction, but no such reversals along 



On-line catalogue at 



http : //www . atnf . csiro . au/research/pulsar/psrcat 
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a north-south axis. This situation must be reasonably localized, because no such reversal is 
seen in the RMs of pulsa rs or of extragalact ic sources at the same Galactic longitude but at 
lower Galactic latitudes (IBrown et al.l 120071). While such "magnetic anomalies" have been 



seen previously in RM (e.g., IBrown fc Taylor! l200l[ ). they tend to be associated with the 



complicated distribution of ionized gas and magnetic fields along the Galactic plane. At 
this relatively high Galactic latitude, with no other extended sources or features anywhere 
in the vicinity, it is hard to see what could cause the required foreground RM distribution. 
It would also need to be a convenient coincidence that the RM distribution of an unrelated 
foreground feature has a symmetry axis that aligns with that defined by the SNR's morphol- 
ogy. We conclude that the RM pattern seen against G296. 5+10.0 is most likely intrinsic to 
the source, and must therefore result from material swept by the expansion of the SNR into 
its surroundings. 



4.3. The RM of a Swept-Up Uniform Ambient Medium 



Kothes fc Brown! (120091 ) have pointed out that ambient gas swept up into a shell by 
an expanding SNR can produce an observable RM signal, and that the resulting systematic 
variations in RM as a function of position around the SNR can be used to infer the in situ 
magnetic field geometry of the gas into which the SNR has exploded. However, there are 
two difficulties arising if the RMs seen toward G296. 5+10.0 are interpreted in this way. 

First, the observed pattern of RMs is completely contrary to expectations. If a SNR ex- 
pands into a well-ordered ambient field with a signi ficant component in the plane of the sky, a 
barre l morphology should result as is observed here ( van der Laanfl~962b ; Whiteoak fc Gardner 
19681 ). but there will be identical RM gradients along each SN R limb, with the gradient 
aligned with the SNR's symmetry axis (IKothes &: Brown! 120091 ) . We observe the exact op- 
posite here: there is a minimal RM gradient along the SNR's symmetry axis, but a large 
gradient across it (Figure [3]). This field geometry is difficult to explain by a model in which 
the SNR sweeps up a uniform magnetic field. 

Second, the magnitude of the RM through the shell implied by this scenario is much 
larger than observed. To illustrate this, we assume that the RM pattern of the SNR is due to 
a geometric effect that is occurring equally in both limbs (deferring a specific interpretation 
of this pattern until §4.41 below). We then compare the RMs of the two limbs to separate 
the foreground RM contribution from that due to material in the SNR shell: the average 
RM between the two limbs is approximately +5 to +10 rad m -2 (which is roughly consistent 
with the foreground contribution estimated in §4.21 above), while the difference between the 
limbs is ~ 40 rad m~ 2 . Thus we estimate that the observed RM contribution of each limb is 
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~ 20 rad m 2 , with the average field directed in opposite directions on each side to produce 
the observed RM pattern. 

Assuming that the synchrotron-emitting and Faraday-rotating regions in the swept-up 
shell are mixed together, the region will not be Faraday thin, but rather will emit a continuous 
spread of superimposed RMs, ranging from a foreground-only RM from material at the edge 
of the SNR closest to Earth to a maximal RM value experienced by polarized signals that 
propagate from the back of the SNR, through the swept-up shell, and then through the 
foreground ISM. The observed RMs and the variation between them are both much smaller 
than the width of the RM transfer function in our observations of ~ 500 rad m -2 . We thus 
cannot resolve the source as Faraday thic k, but instead will measure an average RM equa l 
to half the full value through the SNR ^Gardner fc Whiteoakl Il96~6l : ISokoloff et aD Il998h . 
Thus in relating the observations to the electron density and magnetic field strength in the 
swept-up limbs, the observed RM of ~ 20 rad m -2 implies an RM through the SNR shell of 
RMism ~ 40 rad m~ 2 . 

To compare this to the expected value of RM^m, we assume that a SNR of radius 
R pc has expanded into a medium of uniform hydrogen density no cm -3 . If the SNR shock 
compresses the material that it sweeps up by a factor X, then the thermal electron density 
within the material shocked by the SNR (assuming 100% ionization of ambient gas by the 
SNR shock) is n\ = Xuq where the subscript 1 denotes the post-shock state. Conservation 
of mass requires that the thickness of the swept-up shell is AR = R/3X, while the maximum 
depth of a sightline through the SNR limb is L « 2{2RAR) 1/2 (assuming AR < R). The 
RM through the entire limb is then 0.81ni_BiL cos 9, where B\ is the magnetic field strength 
in the post-shock gas, and 6 is the angle of the post-shock field to the line of sight. The 
maximal RM for the swept-up material, for 6 = 0, is then RM^m ~ 1.3X 1 ' 2 noBiR rad m~ 2 . 

At a distance of 2.1 kpc, a representative radius for the SNR (measured from the geo- 
metric center to the brightest regions of the limbs to the east and west) i s R ~ 20 pc. H I an d 
X-ray observations of the SNR imply no ~ 0.1 cm -3 (IDubner et al.lll986l ; lMatsui et al.lll988l ). 



while equipartition bet ween relativistic particles and magnetic fields in the em itting region 
(e.g., Pacholczyklll970l ) implies strength in the SNR shell of B\ 25 — 50 «G. iRoger et al. 



( 119881 ) argue that SNR G296. 5+10.0 is in the adiabatic (Sedov- Taylor) phase of evolution, 
for which X = 4; the implied shell thickness AR = R/12 is consistent with the observed SNR 
morphology. We then predict RMism ~ 130 — 260 rad m~ 2 , which is substantially larger than 
the value RM^m ~ 40 rad m~ 2 required by observations. From the above considerations, we 
conclude that an ambient uniform medium cannot explain either the strength or pattern of 
the observed RMs. 
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4.4. The RM of the Swept-Up Progenitor Wind 



An alternative explanation for the RMs observed in G296. 5+10.0 is Faraday rotation 
through the magneto-ionic medium of the swept-up stellar wind from the supernova progen- 
itor. At large distances from the stellar surface, stellar winds are expected to have largely 
toroidal fields; if such a field geometry is preserved after material has been swept up by the 
supernova explosion, this could potentially produce an RM pattern that was negative on one 
limb of the SNR and positive on the other. We now consider this possibility in detail. 



The presence of a central neutron star in G296. 5+10.0 (e.g., iHelfand fc Becker! Il984j ; 
Mereghetti et al.lll996[ ) demonstrates that the SNR resulted from the core-collapse of a mas- 
sive star. These stars have powerful winds that can result in substantial mass loss well before 
the supernova explosion, and wh ich can extend for many parsecs from the star (ILozinskaya 
19921 ; iGarcia-Segura et al.lll996aU bl) It is thus both theoretically expected and well-established 
from observations that many SNRs are propagating through their progenitor stel lar winds, 
with consequent implications for an SNRs dynamics and emiss ion properties ( IChevalier 



19821 ; IChevalier fc Liangj|l989l ; iDwarkadad 120071 ; IWeiler et al.H????f). Indeed , such a possibil- 



ity has been specifically suggested for G296. 5+10.0 by IStorey et al.l (119921 ). to explain the 
high degree of mirror symmetry between the two limbs. 



It is a 
fie lds (e.g. 



so well-esta blished that at least some massive stars harbor strong magnetic 



(I Weber &: Davis 



1967 



Wadd 120031). and that thes e stars correspondingly la unch magnetized winds 



Ignace et al.lll998l : lud-Doula fc Owockill2002h. The effects on an SNR 



expanding into a magnetized wind were discussed by IChevalier &: Luol (119941 ). but the result- 
ing RM of the swept-up wind was not addressed. Here we provide a first simple calculation 
of the likely Faraday rotation signature of this interaction; a subsequent study will further 
explore this phenomenon through a detailed magnetohydrodynamical calculation. 

As noted above, an azimuthal magnetic field naturally produces Faraday rotation that 
is positive on one side of the center and negative on the other. For this field to have been 
produced by the stellar wind of the progenitor, we need to ass ume that the symm etry axis 
of the SNR corresponds to the spin-axis of the progenitor (see lStorey et al.lll992l ) and that 
the field pattern in the wind is preserved when it is swept up by the SNR. 

However, we note that the simplest geometry for a magnetized stellar wind is a split 
monopole, in which the magnetic field configurations in the n orthern and southern hem i- 
spheres have opposite polarities, separated by a current sheet ( ]Sakurailll985l ; iMestell 120031 ) . 
For such a wind, the azimuthal component of the pre-supernova magnetic field will be di- 
rected in opposite directions above and below the equatorial plane of the progenitor. The 
simplistic expectation is that a swept-up shell of such material would show a quadrupolar 
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pattern of RMs with, for example, positive RMs to the northwest and southeast, and negative 
RMs to the northeast and southwest. This is not the pattern that we observe in Figure [TJ 

For a magnetized swept-up wind to be a viable explanation of the observed RMs, we 
thus require a higher-order field in the progenitor wind, such a s an aligned quadrupole. 



Such a wind, as may occur in at least some supernova progenitors ( jThompson fc Landstreet 



19851 ). can have an azimuthal field of one polarity in a broad equatorial region covering both 
northern and southern sides of the equatorial plane, with regions of alternating polarity 
toward the poles. Since the barrel morphology of SNR G296. 5+10.0 is dominated by two 
limbs on either side of the symmetry axis, this is a geometry that can potentially match the 
observed RM pattern. Further investigations of this wind geometry and its interaction with 
a SNR shock are needed, but for the purposes of the discussion below, we assume that this 
or other simple wind patterns can produce the RMs seen in our data, and we now consider 
whether the magnitudes of the RMs are consistent with expectations from a swept-up stellar 
wind. 

For a spherical stellar wind of mass-loss rate M and asymptotic velocity Voq, conservation 
of mass implies that the number density n{r) of gas at a radius r from the center of the star 
(where r lies outside the wind acceleration region) is: 

M 

< r ) = 7-^T7 ' 2 

Airr z V 00 mH 

where m# is the mass of a hydrogen atom. 

For a rotating star with a magnetic field, the (presumed ionized) stellar wind is channeled 
by the magnetic field lines i.e. the plasma is frozen into the field. Magnetic field lines leave 
the surface of the star radially and rapidly wind up into a spiral shape at radii greater than 
the maximum radius of co-rotation. For a star of radius with an equatorial rotational 
velocity V rot and assuming symmetry above and below the equatorial plane, the strength, 



B^jr ), of the azimuthal component of the stellar wind's magnetic field (ILamers &: Cassinelli 



1999|) is: 

for r ^> i?*, where -B* is the magnetic field strength at r = i?*. 

We now consider the result of an SNR expanding through this wind, so that the gas 
and magnetic field distributions described by these equations are swept up into a thin shell 
of radius R and thickness AR. For the purposes of the present discussion, we make the 
simplifying assumption that after this material is swept up, the shocked material within the 
shell has a uniform density n\ and a uniform magnetic field strength Bi, as for the simpler 
case discussed in §4.31 above. 
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The total number of atoms in the wind, N to t, swept up by the SNR is: 



R 2 , MR 



Ntot = / n(r)A7rr 2 dr = — (4) 

for r > R*. 

Conservation of mass requires that N tot = 47iR 2 ARni, so that: 

ni = AnR^m, (^) = ^ (^) " (5) 

We can similarly consider the total magnetic flux, $ to4 , passing through a meridional (con- 
stant $) area element that extends in radius from i?* to R and spans Aip in colatitude: 



®tot= / B^rdr = B*R*RAiP ' 
Jr* 



rot 



oo 



(6) 



where ip is the polar angle of the spherical co-ordinate system. 



As the material is ejected radially, the meridional area is radially compressed but the polar 
angle is unchanged before and after the compression. 

Assuming conservation of flux after the SNR sweeps up this material, we require = 
BiRARAip and hence: 

R* \ I Vrot \ n / t->\ / R 



We are now in a position to estimate the RM produced by the swept-up wind, 

RM^ = 0.81 /(£ygM)g) cos eradm - 2 (8) 

where the integral is along a path-length dl through the shell, 9 is the angle of B\ to the line 
of sight, and where we have assumed that the swept-up wind is fully ionized. As in §4.3[ the 
maximum path length through the limb is L « 2{2RAR) l l 2 . Through the limbs of the SNR, 
a swept-up azimuthal field will be oriented such that 9 ps 0. Substituting Equations ([5]) and 
((7j) into Equation (|H|), and writing x = Kot/V'oo, we derive: 



RM wind = 3.6 x 1(T 



R 



3/2 



ar) 

■r \ I m \ ( n v 2 ( v„ 



0.1/ V 1O- 6 M yr- 1 / V20 pc J V1000 km s 



£0 


1 / R* 




' V1OO# 




-i 


-i) 


rad m" 



(9) 
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This equation clearly implies that a wind from a main-sequence or Wolf-Rayet supernova 
progenitor will not produce significant Faraday rotation when it is swept up, because the 
stellar radius is small (< lOi? ), the mass- loss ra te is low (M < 10~ 5 — 10~ 6 M M yr" 1 ) and 



the wind veloci ty is high (V^ > 1000 km s x ) (see I Chevalier &: Liangill989l ; lDwarkadasll2007 



Crowtherll20071 . and references therein). However, if we adopt representative values for a red 
supergiant and it s wind, B* = 500 G R* = 300i^ x = 0.1, M = 1 x 10~ 5 M^ yr" 1 and 
= 30 km s" 1 jGrav k Tonerl ll987l: lDorscdl2004[lviemmings et alil2005l : Ivan Belle et al 



20091 : Ivan Loonl 120091 ). we predict: 



RM 



wind 



u 



—) 

ARJ 



3/2 



rad m 



(10) 



For a shell thickness R/AR fh 8 as indicated by the radio morphology of the SNR, this yields 
RM win d ~ 40 rad m -2 , which is a factor of two larger than the observed RM, as required (see 
discussion about Faraday depth effects in ^4.3[) . We note that a red supergiant wind flowing 
outward for ~ 1 Myr (e.g. IWoosley et al.ll2002l ) at Voo fx 30 km s _1 will extend ~ 30 pc from 
the star, so that an SNR with a current radius R ~ 20 pc would still be propagating through 
this wind. 

We caution that Equation ([9]) has many free parameters, and that supernova progenitors 
show a very wide range of properties. Thus we cannot conclusively demonstrate or rule out 
that the RM seen against the SNR matches the expected Faraday rotation from its progenitor 
wind. However, having ruled out other reasonable alternatives for the observed RM pattern 
in § §4.21 & I4.3[ and having shown that the RM for at least one set of reasonable wind 
parameters can match the observations, we conclude that an azimuthal magnetic field in 
the stellar wind of a red supergiant progenitor can potentially explain both the sign and 
magnitude of the RMs seen against SNR G296. 5+10.0. 

As an additional test for this scenario of a swept-up red supergiant wind, we now consider 
whether the barrel-shape of G296. 5+10.0 is consistent with our proposed model. In order to 
do this, we will examine the two possible mechanisms by which a magnetized red supergiant 
progenitor star can produce a bilateral SNR. These are: (i) expansion of the supernova into 
a magnetized stellar wind and (ii) expansio n of the supernov a into a stellar wind with a 



biconical outflow aligned with the spin axis ( Manchester 



1983). 



In the first case, the circumstel lar material is distributed axisymmetrically by the action 
of the magnetic field (jSakurailll985l ) and is swept up by the expanding supernova shock into 
a barrel-shaped SNR. Similarly, in the second case the supernova sweeps up material in all 
directions except for a bipolar outflow cavity, which has been largely evacuated of material 
due to the high wind velocity along this axis. Under either of these schemes, the result is a 
remnant with two bright limbs (one at either side of the spin axis of the star) and two null 
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points along the axis. 



These two models for the creation of barrel-shaped SNRs both involve a supernova 
expanding through an axisymmetric distribution of circumstellar material. A prominent 
example of such a system is the young barrel-shaped supernova remnant of SN1987A. This 
remnant is understood to have been shaped by the interaction of the supernova shock with 
the anisotropic circumstellar me dium caused by the red and blue supergiant phases of the 
progenitor star ( ]Dwarkadasll2007l ) . We argue that a supernova expanding into a magnetically- 
driven stellar wind is the most plausible sce nario and explains t he structure of G296. 5+10.0. 
This notion is supported by the findings of IStorey et al.l ( 119921 ). who argued that the small- 
scale correlations in structure on either side of the SNR originate from changes in the outflow 
from the progenitor star rather than an external ISM. 

An add i tional consideration is the large axial ratio of SNR G296. 5+10.0. Although 
Roger et al.l (119881 ) attributed the elongation of the limbs to coupling of material with the 
magnetic field of the ambient ISM, a model involving a SNR expanding into a magetized 
stellar wind can just as readily explain this phenomenon. The distribution of circumstellar 
material resulting from a magnetized progenitor stellar wind will result in a large fraction 
of material being deposited in the magnetic equatorial plane and progressively less mass 
distributed towards the polar regions. A spherically symmetric supernova explosion occuring 
within such a medium will result in a remnant with some 'straightening' of the limbs, because 
when the shell is being swept up, the pressure experienced by the SNR shell encountering 
the equatorial CSM will be greater than the pressure encountered by that part encountering 
less dense material towards the polar regions (IBlondin et al.lll996 ). Th is picture is consistent 
with the small-scale structure symmetry noted by IStorey et al.l (119921 ) in G296. 5+10.0. 



A final piece of th e puzzle is the radio polarization observation of SNR G296. 5+10.0 by 
Milne fc Hayned ( 11994] ) . By observing at higher radio frequencies than we present here, they 
were able to infer the intrinsic orientation of the polarization vectors in the SNR's synchrotron 
radiation, before this emission is Faraday rotated in foreground thermal gas. They found 
that the magnetic field orientation of the synchrotron-emitting regions has a projection on 
the sky which is largely tangential to the rim of the SNR. This can be explained in the 
context of a magnetized progenitor wind if we consider a progenitor star with its spin axis 
tilted with respect to the plane of the sky and with an azimuthal magnetic field. Such 
an alignment would produce a supernova remnant with linearly polarized emission that is 
aligned with the sky-plane magnetic field direction tangential to the remnant limbs, as well 
as a line-of-sight magnetic field that reversed direction between the two separate limbs, as 
observed from our Faraday rotation measurements. 
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5. Conclusions 

Using rotation measure synthesis, we have made Faraday rotation maps of the supernova 
remnant G296. 5+10.0. The limbs of the remnant were found to have RMs of opposite 
signs and separated by approximately 40 rad m -2 , implying that the magnetic field in the 
foreground to the polarized SNR emission is oriented in opposite directions on either side 
of the SNR's symmetry axis. By comparing the published RMs of pulsars observed through 
similar sight-lines through the Galaxy, we demonstrated that the spread of RMs within the 
remnant is inconsistent with Faraday rotation by the foreground magneto-ionized medium. 
We also considered whether a swept-up ambient interstellar medium can explain the RM 
distribution in G296. 5+10.0, but found a large discrepancy between the observed rotation 
measures and those predicted by theory. Finally, we found that the swept-up wind of a red 
supergiant progenitor star could explain the observed magnitude and distribution of RMs in 
the remnant, as well as the linear polarization found by other authors. This would require 
an aligned quadropole or similar magnetic field morpology in the wind of the progenitor star. 

Future observational studies should use RM synthesis to map Faraday rotation in a larger 
sample of Galactic supernova remnants. Focusing on remnants that lie at relatively large 
Galactic latitudes will minimize the effect of the Galactic foreground and enable an accurate 
determination of the intrinsic properties of the magnetized plasma within the remnants. If 
similar morphologies are found, we will be able to gain a deeper understanding about the 
mass-loss histories of progenitor stars prior to their supernova events. If the Faraday rotation 
in the remna nt represents the signa ture of a magnetized stellar wind, then the method 



suggested by iKothes fc Brown! (120091 ) for studying ambient magnetic fields in the Galaxy 
using SNRs is not only a powerful in situ probe of magnetic fields, but also of SNR progenitor 
winds. By studying a larger sample of sources we can hope to understand the physical origin 
of magnetic fields swept up in supernova remnants. In parallel, magnetohydrodynamical 
simulations of red supergiant winds should also be employed to test the validity of our 
model of a swept-up magnetized stellar wind in G296. 5+10.0. 
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Fig. 1 — ATCA images of SNR G296.5+10.0 at 1.4 GHz. (a) Total intensity, shown with a 
greyscale ranging between —2 and +100 mJy beam -1 . The blue ellipse at the lower left is the 
angular resolution of the image, 3'.3 x 1'.8, at a position angle of —142° north through east; 
(b) Linearly polarized intensity, with greyscale ranging between 0.5 and 70 mJy beam" 1 ; (c) 
Rotation measure plotted every five pixels (150"), with filled boxes representing positive RMs 
and unfilled boxes representing negative RMs. The linear dimensions of each box indicate the 
magnitude of the rotation measure. The range of box sizes corresponds to — RM — ranging 
from to 128 rad m~ 2 . Contours are the Stokes-/ image at +20 and +80 mJy beam _1 ;(d) 
An image of RM, with the color scale indicating RMs between —50 rad m -2 (red) and +70 
rad m~ 2 (blue) as indicated to the right of the image. 
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Fig. 2. — The RM synthesis spectra of representative pixels in the eastern (blue) and western 
(red) limbs of SNR G296. 5+10.0. Spectral deconvolution has been applied to remove the 
sidelobes caused by incomplete sampling in A 2 space. 
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Fig. 3. — Plots of RM against right ascension (top) and against declination (bottom) in 
G296. 5+10.0 for the RM distribution in Figure 1 (c); red points are in the western limb and 
blue points are in the eastern limb. The mean RM for each limb is displayed as a dashed 
line. 



